Effects of Inhomogeneity on the Spectrum of the Mott-Insulator State 
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We investigate the existence of quantum quasi phase transitions for an ensemble of uhracold 
bosons in a one-dimensional optical lattice, performing exact diagonalizations of the Bose-Hubbard 
Hamiltonian. When an external parabolic potential is added to the system quasi phase transitions 
are induced by the competition of on-site mean-field energy, hopping energy, and energy offset among 
lattice sites due to the external potential and lead to the coexistence of regions of particle localization 
and delocalization in the lattice. We clarify the microscopic mechanisms responsible for these quasi 
phase transitions as a function of the depth of the external potential when the on-site mean-field 
energy is large compared to the hopping energy. In particular, we show that a model Hamiltonian 
involving a few Fock states can describe the behavior of energy gap, mean particle numbers per 
site, and number fluctuations per site almost quantitatively. The role of symmetry on the gap 
as a function of the depth of the external trapping potential is elucidated. We discuss possible 
experimental signatures of quasi phase transitions studying the single particle density matrix and 
explain microscopically the occurrence of local maxima in the momentum distribution. The role of 
a thermal population of the excited states on the momentum distribution is discussed. 

PACS numbers: 03.75.Kk, 03.75.Lm, 05.30.Jp, 73.43.Nq 



I. INTRODUCTION 



The theory of quantum phase transitions for bosonic 
atoms in an optical lattice has received increased interest 
due to the recent observation of the Mott-Insulator tran- 
sition created by loading a Bose-Einstein condensate . 
Optical lattices are created by pairs of counterpropagat- 
ing laser beams, which lead to spatially dependent light 
shifts for the atoms. Early theoretical work showed the 
presence of quantum phase transitions resulting from the 
interplay of on-site interactions and hopping of atoms 
between lattice sites which can be controlled by vary- 
ing the intensity of the lattice laser beams Numer- 
ical Monte Carlo techniques have been succesfuUy used 
to locate the superfluid to Mott-Insulator transition for 
commensurate number of particles and wells The 
addition of an external harmonic potential leads to the 
coexistence of domains where particles are localized and 
domains where particles are delocalized 0, |3 . Reference 
shows that the distinction between commensurate and 
incommensurate filling of the lattice disappears in pres- 
ence of the external trap. Suggestions for the detection 
of transitions in a trap have been proposed . An exact 
numerical study for a limited number of atoms and wells 
shed light on the effect of patterned potentials 0- The 
influence of the excited-state spectrum on superfluidity 
is studied in Refs. 

In this paper we study atoms in tightly confining one- 
dimensional lattices in presence of a quadratic external 
potential. Such one-dimensional lattices can be real- 
ized by independent control of laser intensities and wave- 
lengths along the three spatial directions. The quadratic 
potential is either generated using a magnetic trap or a 
patterned optical lattice, where the longest-wavelength 
lattice can be approximated as being harmonic (IQi] . Ex- 



act numerical results, based on the second-quantized 
single-mode Bose-Hubbard (BH) Hamiltonian jj| , are ob- 
tained for a limited number of particles and wells. The 
paper is organized as follows: first we review the spec- 
tral properties of the homogeneous system and discuss 
the character of the ground-state wave function. Then 
we clarify the microscopic mechanisms responsible for 
external-trap-induced phase transitions when the on-site 
mean-field energy is large compared to the hopping en- 
ergy. We refer to these transitions as quasi phase tran- 
sitions both because our study involves finite number of 
particles and wells and our study investigates the simul- 
taneous presence of localized and delocalized domains 
caused by the presence of the external trapping poten- 
tial. The behavior of observables like energy gap, mean 
particle numbers and number fluctuations in the differ- 
ent phases is explained using a model Hamiltonian. In 
the strong mean-field regime the calculations are read- 
ily generalized to large commensurate numbers of par- 
ticles and wells. Finally, we discuss possible experimen- 
tal verifications of the transitions in current experiments, 
the effects of thermal population of excited states on the 
single-particle momentum distribution, and conclude. 
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II. THE BOSE-HUBBARD HAMILTONIAN 
A. Bose-Hubbard Hamiltonian 

The one-dimensional Bose-Hubbard Hamiltonian is 
given by 



^ 1 1 t 



(1) 



where the integer i labels the sites, and a\ are the 
annihilation and creation operators at site i respectively, 
J and U are the site-independent hopping constant and 
mean-field energy, respectively, and the are energy off- 
sets that depend on index i and describe the external 
trapping potential. Particle hopping is limited to near- 
est neighbors. The number of wells and particles is M 
and N , respectively. For this paper we assume that M 
is odd and A'' is an integer multiple of M . The har- 
monic external trap is then given by — ed{i/ L)"^ , where 
L — [M — l)/2 is integer and i runs from —L to L. In 
this way = cl ~ to is the "trap depth" between the 
outer-most, i — ±L, and central, i = 0, well. The Hamil- 
tonian is invariant under spatial reflection / around the 
central site z = 0. 



B. Basis set 

We solve the many-body Hamiltonian in the occupa- 
tion number represention or Fock-state basis. That is, wc 
define a symmetrized basis function by \{q-L ■ ■ ■ Ql}^) = 
{\q-L ■ ■ ■ 9l)±|'?l ■ ■ • 9-l))/-\/2, where qt is the number of 
atoms in site i. A symmetrized basis function is an eigen- 
state of the reflection operator I, with eigenvalue ±1. 
The number of basis states is (Af + - 1)!/((M- 1)!7V!), 
where the exclamation mark denotes the factorial. In 
this basis the diagonal matrix elements of the Hamil- 
tonian are given by J^i Uqi{qi — l)/2 + e^gi. The non- 
zero off-diagonal matrix elements are solely determined 
by the hopping interaction. The eigensystem is solved 
with standard linear algebra techniques. Eigenstates are 
classifled by reflection symmetry and are named "ger- 
ade" and "ungerade" if they are even or odd under the 
symmetry, respectively. 



III. HOMOGENEOUS LATTICE 



bands are visible. This regime is known as Mott-regime. 
The non-degenerate ground state wave function is to a 
good approximation given by the symmetric Fock state 
|{111111111}+) = llllllllll). Eigenstates belonging to 
the second band are principally dominated by the 72 Fock 
states with site occupation 0,2,1,1,1,1,1,1, and 1. These 
Fock states have the same diagonal matrix element. The 
third band, of which we show only the eigenstate of low- 
est energy, consists of states whose contribution comes 
from 0,2, 0,2, 1,1, 1,1, and 1. The energy difference between 
successive bands is therefore approximately U , the mean- 
field energy needed for a particle to hop to an already oc- 
cupied site. For U / J < b no bands are discernable, cor- 
responding to the superfluid regime. No single Fock state 
dominates the ground and excited-state wave functions. 
Reference Q has shown that the transition between the 
superfluid and the Mott-state occurs at {U/J)c — 4.65, 
for N/M = 1. 



B. The Mott-Insulator state 

We investigated the nature of the ground state in the 
Mott-insulator regime in more detail. For C// J 3> {U/ J)c 
the dominant conflguration is |111111111). Particle hop- 
ping mixes so-called excitonic states into the ground state 
wave function, 111. Here excitonic states are Fock states 
in which one site has no particles, one of its two near- 
est neighbors has two particles, and other sites have sin- 
gle particle occupation. Particle-hole states in which the 
doubly occupied and the empty site are not nearest neigh- 
bors do not couple to the |111111111) state. In second 
order perturbation theory the excitonic mixing lowers 
the ground-state energy by 4(M — /U . There are 
2(M — 1) excitonic states. 

Figure|21shows the binding energy for the exact ground 
state and for the second order model. Note that the ex- 
pectation value of the Hamiltonian over the dominant 
Fock-state |111111111) is zero. The two curves agree 
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A. Energy spectrum 




Figure n shows the lowest 74 eigenenergies of the spec- 
trum for a homogeneous lattice, = 0, and commensu- 
rate filling, N = M = Two physically distinguishable 
regimes can be discerned. For U/ J > 5 three energy 



FIG. 1: Spectrum of the BH Hamiltonian in a homogeneous 
lattice and N = M = 9 as a function of U/J. The first two 
bands and the lowest state of the third band of the complete 
spectrum are shown. The first band has a single state. 
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for U/J>b. In other words the ground state is well 
described by perturbation theory even close to the tran- 
sition point. 

The zero-temperature single-particle density matrix is 
defined as pij = {a\aj), where the average is over the 
ground-state wave function. Mixing of excitonic states 
into the ground state wave function gives rise to non- 
zero off-diagonal elements of ptj. For \J j J {U/J)c 
non-zero elements occur for nearest-neighbour sites and 
are on the order of J/U. The diagonal density-matrix 
elements pu = (rii) , where the operator rii — a|ai, de- 
fine the mean number of particles at site i. The quantity 
CTj = y/ {nf) — (rii)^ defines the number fiuctuations at 
site i. Notice that (nf) > {ni)^. In a homogeneous lat- 
tice (rii) = N/M for all U and J. For the parameters 
in Fig. n the number fluctuations are of order unity for 
U/J ^ {U/J)c and are suppressed for U/J ^ {U/J)c- 
Gi <^ {rii) implies particle localization. We discuss rela- 
tionships between the single-particle density matrix and 
experimental observables toward the end of this paper. 



IV. INHOMOGENEOUS LATTICE 

For the remainder of this paper we investigate the 
regime U/J » {U/ J)c by adding an additional parabolic 
trap. Zero-temperature quasi phase transitions are in- 
duced by varying the depth of the parabolic trap and 
correspond to a redistribution of particles in the lattice. 
A microscopic explanation of these quasi phase transi- 
tions in terms of a few dominant Fock states is given. 
These results will largely be done by setting U/J = 30 
(see dotted line in Fig.QJ and varying e^. 



A. Observations 

Figuresl^land ^show the mean number of particles and 
number fluctuations of the ground state as a function of 
trap depth td for N = M = Q and U/ J = 30, respectively. 
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FIG. 2: Binding energy of the ground state for the exact cal- 
culation (full line) and the perturbative model (dashed line) 
as a function of U/J, ior M = N = 9. 
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FIG. 3: Mean particle number (rn) at site i for 9 particles, 9 
sites, and U/J = 30 as a function of trap depth e^/J. Curves 
are labeled by their site index i. Regions of particle local- 
ization are characterized by integer occupation, while non- 
integer occupation signals particle derealization. 

The sum of {rii) over all sites equals N = 9. Rapid 
changes or quasi phase transitions are observed in both 
{rii) and (7^. For < U, {rii) is one for all sites on the 
scale of the Figure, while the number fluctuations are 
about 0.1 particles and are slowly increasing. 

Near Cd — U, jumps in {rii) and ai are apparent. The 
mean particle number in the outer-most wells drops to 
0.5 and then to zero, while simultaneously {rii) increases 
for the central and ±1 wells. In fact, the increase oc- 
curs first for the central well and then for the ±1 wells. 
The number fiuctuations for the outer-most wells jump 
to about 0.5 particles and then fall back to 0.1 particles. 
Fluctuations for the ±1 wells increase to about 0.4 parti- 
cles and then to 0.5 particles, ctq jumps to « 0.5 particles 
and then slowly decreases as a function of ed- 

The observations signal two quasi phase transitions in 
the system. The first one corresponds to a spatial rear- 
rangement of the atoms, which creates domains of par- 
ticle delocalization, with non-integer occupation, at the 
center and at the edges of the lattice, and domains of 
particle localization, with integer occupation, in between. 
After the second quasi phase transition particles are de- 
localized in the three central wells only. 

Jumps in mean number of particles for larger indi- 
cate successive transitions. For example, the next tran- 
sition occurs at about ed — 62 J. Redistribution from 
sites \i\ < L towards the central wells occurs. The first 
two transitions are of particular importance as they oc- 
cur for any N and M, with N/M integer, and contain 
all necessary physics that is needed to understand parti- 
cle rearrangements for other quasi phase transitions. We 
therefore discuss them in detail. 



B. Energy spectrum 

Figure [S] shows the spectrum of the Hamiltonian for 
N = M = 9 and U/J = 30 as a function of the trap 
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depth. The expectation value of the operator cirii 
over one of the symmetrized Fock states |{011121111}^) 
has been subtracted to facihtate the interpretation of the 
data. Energies of eigenstates are shown for states that 
originate from the lowest four bands in the absence of 
the external trap. Energy differences among the bands 
become smaller for increasing e^. The graph has series of 
level crossings among energies of all four bands. Cross- 
ings with levels of the fourth band occur for Cd > 25 J, 
and play a minor role in the first two quasi phase transi- 
tions. Many of the crossings are narrow avoided crossings 
between eigenstates of the same symmetry under reflec- 
tion through the center of the lattice. For ej, — U the 
slope of the ground state energy changes abruptly. In 
fact rapid changes occur at ~ 28 J and ~ 32 J. They 
correspond to avoided crossings with the lowest-energy 
gerade state of the second and third band, respectively. 
The lowest curve of the second band contains two nearly 
degenerate states which are indistinguishable on the scale 
of the graph. One of these states is gerade and the other 
one is ungerade. The two avoided crossings lead to the 
first two quasi phase transitions discussed in reference to 

Figs. El and H 

The inset of Fig. O shows the energy gap Egap as a 
function of the trap depth, where the energy gap is the 
difference between the energy of the first-excited and of 
the ground eigenstates. For < U the energy gap de- 
creases with increasing e^. After the first quasi phase 
transition Egap is very small compared to the energy scale 
J. For ed > 32 J the gap is of order J. The ground and 
the first-excited states are gerade and ungerade for all ed, 
respectively. 



C. Effective model 

The behavior of the ground-state energy, energy gap, 
and of (rii) and ai as a function of Cd up to Cd < 60 J 
can be explained in terms of an effective model involv- 
ing six symmetrized Fock states. Table lists the six 
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FIG. 4: Number fluctuations ai at site i for 9 particles, 9 sites, 
and U/ J — 30 as & function of trap depth ed/J. Fluctuations 
of order 1/2 signal the sharing of a particle between two sites. 



states, together with their zeroth-order energies obtained 
by evaluating the diagonal matrix elements of the origi- 
nal BH Hamiltonian for these six states. They belong to 
the lowest three bands for ed = 0. Within each band, the 
states in Table |2 have the lowest value for the diagonal 
matrix element of the BH Hamiltonian. Exact numerical 
calculations confirm that these Fock states dominate the 
expansion of the eigenstates near the quasi phase transi- 
tions at ed — U. 

In the effective model, ignoring for a moment off- 
diagonal couplings and setting ed — 0, |1111111111) 
and |{011121111}^) are the ground and first two excited 
states, respectively. For > their energies depend 
linearly on the trap depth and cross at e^j — U . For 
ed > U the |{011121111}^) states are the new ground 
states and a quasi phase transition has occurred. It is 
energetically favorable for the system to move an atom 
from the outer to the central well. At ed = (16/15)C/ 
the states |{011121111}±) cross |{011122110})^. A sec- 
ond quasi phase transition occurs. Exact calculations 
for U — 30 J show that quasi phase transitions occur at 
ed ~ 27. 9J and ed ~ 31. 3J, consistent with the effective 
model. 

There is no direct coupling in the BH Hamiltonian 
between the first five states of Table |3 Coupling be- 
tween these states requires at least a third-order pro- 
cess in the hopping interaction. On the other hand, 
the hopping interaction directly couples the gerade state 
|{011122110})+ and the sixth state |011212110) belong- 
ing to the third band. The matrix element is —2 J. This 
coupling shifts the second quasi phase transition to a 
somewhat smaller trap depth e* . The energy gap within 
the model is U—ed for < U. For U < ed < e* the gap is 
zero, while for e* < ed it is never larger than w 2 J. Exact 
calculations show an almost quantitative agreement for 
Egap despite the simplicity of the model. The non-zero 
gap for U < ed < e* in Fig. [S] is due to weak couphng to 
states, which are not included in the model. 



D. Inversion symmetry and observables 

An energy gap that is large compared to J can occur 
when in the model the ungerade companion of the ger- 
ade state with the smallest diagonal matrix element does 
not exist. An example of this situation is the ed region 
where the ground state is |111111111) and Egap ~ U — ed- 
This requirement is necessary but not sufhcient. In fact 
when in addition to not having an ungerade companion 
the particle occupation of the ground state wavefunction 
is not the same in all sites, the energy gap is determined 
by energy differences with Fock states that belong to the 
same band in the absence of the external trap. The en- 
ergy gap is then proportional to ed- This splitting can be 
large or small compared to J. An example for nine sites 
occurs for the ed > W region where the ground state is 
|001232100) and the first excited state is approximately 
|{001322100}-). The energy gap is then « ed/L^ which 
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TABLE I: Table of the six symmetrized Fock states of the ef- 
fective model Hamiltonian for commensurate filling q = N/M, 
q integer. The model predicts the locations of the first two 
quasi phase transitions as a function of e^. The band index, 
B, for each state a.t ed — and the diagonal matrix elements 
of the effective Hamiltonian are given. Here, g< = q — 1, 
q> = g -I- 1, and F = Ylf^-iL-i) i'^ / ■ The notation (s)i 
implies s particles in site i. Double dots imply sites with g 
particles. 

Energy 



B 


State 


1 




. • ■(<?)-! (9)0(9)1 ■ 




2 


m<)-L 


■■{q)-i{(l>)o{q)i 






m<)-L 


••(g)-i(<7>)o(g)i 




3 


m<)^L- 


(g)-i('7>)o(g>)i 


il<)LV) 




m<)-L- 


('?)-i('?>)o('7>)i 


•(9<)i}") 




\{1<)-L 


■('7>)-i(g)o(g>)i 


■•(g<)^) 



(r+2)6d 

U+{T+l)ed 

2[/+ F+l/L^ ea 
2U+ {v+2/L^) 



is large compared to J. 

The energy gap is much smaller than or on the order 
of J when the ungerade state exists. In this case, the two 
symmetrized Fock states have identical diagonal matrix 
elements. In fact, the energy gap is much smaller than 
J when within the model there is no first-order coupling 
between the two gerade states with the smallest diagonal 
matrix elements. For a gap of order J direct coupling 
exists between the two gerade states. An example of a 
gap that is small compared to J occurs for the trap depth 
range where the ground state is |{011121111}+) and the 
first excited state is |{011121111}^). The gap is of order 
J when the ground state is |{011122110})+. 

The observables {rii) and a±i in Figs. El and 01 de- 
scribe ground-state properties of the system, and their 
behavior through quasi phase transitions can be qual- 
itatively explained from the gerade Fock states of Ta- 
ble Ql Integer (n^) for all i ( and therefore {ni) — {n^i) 
) and suppressed Ui corresponds to a situation in which 
the dominant Fock state of the ground state is of the 
form llllllllll) or |001232100), where the correspond- 
ing ungerade state does not exist. Non-integer {rii) and 
Gi of order (rii) in two or more sites signals particle de- 
localization between those sites. Examples of delocaliza- 
tion occur for [/ < < e* in sites ±L,±1,0 and for 
td > e* in sites ±1,0. The approximate ground state is 
|{011121111}+) and |{011122110})+, respectively. 

Particle delocalization can be classified in terms of the 
distance between the sites in which it occurs. If the do- 
mains of delocalization are not connected, that is, the 
sites are not nearest neighbors, {rii) and Ui assume ap- 
proximately half-integer values. Half-integer values are a 
consequence of the inversion symmetry of the problem. 
Disconnected delocalization is observed between sites ±L 
for U < ed < €*■ A "Schrodinger-cat" is created between 
those sites. If the domains where particle are delocalized 
are connected, (rii) and at can assume any value. Con- 
nected domains can only be centered around the bottom 
of the parabolic trap, where the differences between 
for neighboring sites are smallest. For U < ed < e* we 
have not only disconnected domains, but also connected 




FIG. 5: Spectrum of the BH-Hamiltonian as a function of 
ed/J, for N = M = 9 and U/J = 30. The expectation value 
of the operator e^ni on the Fock state |{011121111}^) has 
been subtracted. The inset shows the energy gap Egap as a 
function of e^/ J on a logarithmic scale. 



delocalization among the three central sites. A model 
description which reproduces the observed mean num- 
bers and fluctuations in the central sites must take into 
account mixing with states other than those in Tabled 
The energy gap is much less than J for disconnected do- 
mains and of order J for connected domains of delocal- 
ization. Particle delocalization among central sites can 
be understood within our model for > e*. 



E. Extensions of the model 

The microscopic analysis of the energy gap and other 
discussed observables can be generalized to different sit- 
uations. Successive transitions at larger ed can be calcu- 
lated using similar perturbation theory. For a number of 
particles and wells which is not commensurate, the spec- 
trum is qualitatively different from the commensurate 
filling when the parabolic trap is absent. For example, if 

= M + 1 and J7/ J » 1 the lowest band has M-levels 
separated by energies of order J as opposed to the single 
ground-state level of the commensurate filling. More- 
over, Ref. j3| recognized that for incommensurate fillings 
the ground state is superfluid for any U. Nevertheless, 
when the trap is present a perturbative discussion with a 
limited number of symmetrized Fock states is still valid. 
For each band the Fock state with the lowest diagonal 
matrix element must be taken into account in analogy to 
the discussion for commensurate filling. 



V. MOMENTUM DISTRIBUTION 

In current experiments, both lattice and external po- 
tentials are switched off and the atomic density is mea- 
sured after a certain time of flight. A spatial inter- 
ference pattern is observed, which is directly related 
to the momentum distribution or Fourier transform of 
the zero-temperature single-particle density matrix be- 
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FIG. 6: Momentum distribution n{k) for iV = M = 9, U/J = 
30 as a function of k. The gray, dotted, full and dashed 
curves are for ed ~ 0, 26J, 30J, and 33J respectively. The 
momentum distribution is periodic with period 2n. 

fore the trapping potentials are turned off 0, Q ■ That is 
the measurement yields 

where the sums are over all lattice sites. It is tempting 
to try to map the momentum distribution to connected 
or disconnected delocalization within the ground state. 
For U/J ^ {U/J)c, where the ground state is well ap- 
proximated by a few Fock states, the connection might 
be easy to establish. 

A. Observations and discussion 

Figure IHl shows the momentum distribution for four 
values of the trap depth with N = M — 9 and U/J = 30. 
Two of the ed values are smaller than 27. 9J, where the 
first quasi phase transition occurs, one value lies between 
27. 9J and 31. 3J, and one value is larger than 31. 3J. 
At ed ~ 31. 3J the second quasi phase transition occurs 
(see Figs. Eland 0}. The two curves with the smallest 
ed are slowly varying and nearly indistinguishable. For 
ed — 30 J rapid oscillations are observed on top of a pro- 
nounced peak centered at fc = 0. At ed = 33 J, after the 
second quasi phase transition, the oscillations disappear, 
while the peak ai k — remains. 

For ed < U , the momentum distribution is nearly in- 
dependent of ed- In this region the ground state is to a 
good approximation given by the Fock state |111111111) 
with weak mixing to excitonic Fock states, as discussed 
in reference to Fig. |5] The momentum distribution of 
the state |111111111) is one, such that all momenta are 
equally populated. The remaining weak long-wavelength 
oscillation in n{k) is due to the excitonic states, where 
particle and hole are in neighboring sites. These states 
give rise to nearest-neighbor coherence or off-diagonal 
non-zero elements of pij for j = i -f 1. Consequently, 
the k dependence is cosine-like with period 27r. 




-X).5 -0.25 0.25 0.5 

k / (27t) 

FIG. 7: Momentum distribution n(fc) at finite temperature 
T for iV = Af = 9, U/J = 12, = 10 as a function of k. 
The full, dashed and dotted curves are for fcsT = QJ, O.IJ, 
IJ, respectively. The momentum distribution is periodic with 
period 2n. 

The occurrence of the first quasi phase transition at 
27.9 J is signaled by the appearence of local maxima at 
roughly integer multiples of {2tt)/{M — 1). The ground- 
state wave function is predominantly determined by the 
Fock state |{011121111}+) and has coherence between 
the outer-most sites which share a particle. Accordingly, 
Pij has non-zero off-diagonal elements for i — —L and 
j = +L, and n{k) has a cos(fc(M— 1)) contribution. The 
peak at fc = is due to coherence from residual particle 
fluctuations in all sites as shown in Fig. ^ In particular, 
Fock states with single or double particle-hole pairs con- 
tribute to non-zero off-diagonal elements of the single- 
particle density matrix. Our six-channel model cannot 
reproduce this feature satisfactorily. 

The disappearance of local maxima signals the second 
transition at 31. 3J. The approximate ground state is 
then 1(011122110}'^), which does not have coherence be- 
tween the outer-most sites as both sites are empty. The 
peak at fc = is again due to particle fluctuations in 
all sites. In particular, Fock states with two and three 
particle-hole pairs and states with three particles in a 
site contribute to the non-zero off-diagonal elements of 
Pij. The momentum distribution is not well reproduced 
in our model. In fact, the model predicts cos(2fc) con- 
tributions to n{k), which produce maxima at A; = ±7r 
reflecting the coherence among the sites i = ±1. Fluc- 
tuations in all sites overwhelm this cos(2A:) contribution. 
Such maxima can therefore only be observed for larger 
U/ J ratios, where fluctuations are further suppressed, 
and the agreement with the few-channel model is im- 
proved. 



B. Detection of transition points 

In parameter regimes where a particle is shared or de- 
localized between sites which are not nearest neighbours, 
it is convenient to deflne a coherence length as the 
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distance between these sites. When such delocaUzation 
is present, n(k) can show osciUations at integer muhiples 
of (27r)/Cc- For example for = 30 J shown in Fig. El 
Cc = M — 1. Detection of these oscillations could be 
used to determine transition points experimentally and 
infer spatial density distributions in the lattice. We ob- 
serve clear oscillations when Cc equals the distance be- 
tween the outer-most occupied sites of the lattice. When 
the coherence length is smaller than this separation the 
interpretation of the interference pattern can be compli- 
cated by the presence of residual number fluctuations in 
inner and outer- most sites. Additional information from 
energy-gap measurements can be used to infer spatial 
atomic density distributions in the lattice. 

C. Finite temperature effects 

Finite temperature might be a limiting factor in the de- 
tection of quasi phase transitions. At finite temperature 
T the momentum distribution is the Fourier transform of 
the thermalyzed density matrix 

where the sum is on all the eigenstates I, () is an aver- 
age over eigenstate Z, and kg is the Boltzmann constant. 
Figure [7| shows the momentum distribution for U/J = 12 
and ed/J = 10 as a function of k for different tempera- 
tures. Fast oscillations are observed for kgT = O.IJ, and 
disappear at fcsT = IJ. In this region the approximate 
ground state is |{011121111}+) as the trap depth is such 
that U < Ed < e* and we are in between the two quasi 
phase transitions discussed in reference to Fig. The 
T = curve of Fig. [7| is to be compared to the = 30 J 
curve of Fig.|^ The = peak is now more pronounced, 
consistent with being less deep into the Mott insulator 
regime for smaller U . Nevertheless, fast oscillations can 
be seen for large |fc|. 

The energy gap is much smaller than J for this value 
of €d and U . In fact for U < < e* Egap << J , as 
discussed for Fig. [S] The first excited state is predom- 
inantly given by |{011121111}"). A particle is shared 
between the outer-most sites. In the contribution of the 
first excited state to the thermally averaged n{k) the 
cos(fc(M — 1)) oscillations are out of phase with respect to 
the corresponding ground-state ones. Consequently, for 
temperatures larger than Egap/ks the oscillations disap- 



pear. In addition, thermal averaging leads to broadening 
of the peak around k — 0. For U/J — 30 and €d ~ 30J, 
the energy gap is few times 10~^J, and thus requiring 
very small temperatures in order to observe the oscilla- 
tions. 

VI. CONCLUSIONS 



In summary, we have shown that quasi phase transi- 
tions at zero temperature of an ensemble of bosons in a 
lattice with a confining potential can be microscopically 
understood. The ground state structure of the Mott 
phase in the homogeneous system has been clarified. In 
particular, we studied the role of "excitonic" Fock states 
in the ground-state wave function. When a parabolic 
trap is added, a hierarchy of changes in the nature 
of the ground state, corresponding to the quasi phase 
transitions, lead to coexistence of regions of particle 
localization and delocalization. These changes have been 
explained in terms of a model Hamiltonian involving 
a finite number of symmetrized Fock states. These 
states have been chosen to have the smallest value of 
the diagonal matrix element of the BH-Hamiltonian. 
The behavior of the energy gap, mean particle number 
per site, and number fluctuations per site have been 
characterized near several quasi phase transitions using 
such model Hamiltonian. A discussion of the momentum 
distribution at zero and finite temperature for different 
depths of the confining potential suggests that it might 
be possible to detect the exact points at which transi- 
tions occur and to relate the interference patterns to 
definite coherences/delocalizations of atoms in distant 
sites. In particular, we found that when coherence is 
established between the outermost sites peaks appear 
in the momentum distribution at integer multiples of 
(27r/Cc), where (c is the distance between those sites. 
The presence of these peaks has been explained using the 
model Hamiltonian. If temperatures smaller than the 
energy gap are achieved, detection of analogous peaks 
in the density distribution of matter-wave interference 
patterns can be used to determine experimentally 
transition points. 
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